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Abstract. Intrinsic toroidal rotation in a tokamak can be driven by turbulent momentum 
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evaluate each contribution using gyrokinetic simulations. We find that the relative importance of 
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negligible contributions due to neoclassical poloidal electric field and poloidal gradients of density 
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1. Introduction 


Significant ion toroidal flow has been observed in many tokamaks without any external 
momentum injection [TJ. The origin of this intrinsic rotation is primarily turbulent 
momentum redistribution. Understanding the intrinsic rotation observed in experiments 
is useful because the toroidal flow stabilizes MHD modes m e] and increases the 
confinement time H El El d- In this paper, we investigate the radial turbulent 
momentum transport due to neoclassical flows, which always exist in tokamak plasmas. 
In [8] it is argued that this type of momentum redistribution dominates over other 
momentum flux drives when the size of the turbulent eddies is sufficiently small. 
Comparisons with experiments J9J indicate that this type of momentum flux may explain 
certain features of intrinsic rotation, such as rotation reversals in. 

There are momentum sources and sinks at the boundary of the confinement region, 
and these sources and sinks likely determine the rotation at the separatrix. The radial 
transport of toroidal angular momentum redistributes the momentum within a tokamak, 
and determines the size and sign of the ion toroidal flow given the boundary condition. 
If there is no external momentum injection, the change of the ion rotation is determined 
by the total momentum conservation equation 


d(riimi\ i) 

dt 


= -V ■ 


Pi +Pe±{ I -bb) +p e |jbb 


+ -J x B, 
c 


( 1 ) 


where n*, m;, and V; are the ion density, mass, and fluid velocity, pj= f d 3 vfirriiW is 
the ion pressure tensor, /j is the ion distribution function, b is the unit vector parallel 
to the magnetic held B, p e j_ and p e y are the electron pressure perpendicular and parallel 
to the magnetic held, I is the identity tensor, c is the speed of light, and J is the current 
density. Here, the small electron inertial term and the small electron viscosity terms are 
neglected. Multiplying the toroidal component of equation (JT|) by major radius R and 
taking the hux surface average (...),/, = (1/U 7 ) f dddp (.../B • S/d) and the average over 
several turbulence characteristic times and length scale result in the equation 


I = -UA(v'n) +1 ((J. v^) r , (2) 

where V' — j dddp (B • V0) -1 is the hux surface differential volume, and n = 
m i((f d 3 vfiR(v ■ <£)(v • V'0))^ ) )r is the radial transport of ion toroidal angular 
momentum. Here, 6 and ip are the poloidal and toroidal angle, respectively, and 
is the poloidal magnetic hux, giving the magnetic held B = IS7p + Vp x V^, where 
/ = B V R and B v is the toroidal magnetic held. The transport average (...)r is used 
for averaging over turbulent phenomena that have much shorter time and length scales 
than those in which transport happens. The radial current in the second term on the 
right hand side vanishes due to quasineutrality, so the change of ion toroidal angular 
momentum is determined by radial momentum transport, 


d_ 

dt 



1_ d_ 


(U'H). 


(3) 
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The radial transport of particle, energy and momentum in a tokamak is caused 
dominantly by microturbulence (turbulent transport) rather than by collisional 
processes (classical and neoclassical transport). Microturbulence in a tokamak is 
generated mainly by drift wave instabilities due to radial gradients of temperature and 
density. The radial momentum transport is more complex than particle and energy 
transport because a symmetry of the turbulence in an up-down symmetric tokamak 
makes the radial momentum transport vanish to lowest order in an expansion in p*, 
where p* = p^/a <C 1, pi is the ion Larmor radius and a is the minor radius muni 
(see section 2). The symmetry of the turbulence is broken by flow and flow shear. The 
effect of the flow and its gradient on the momentum transport can be linearized in the 
low flow (sub-sonic flow) regime, giving an advective term and a diffusive term that are 
proportional to the flow and the flow shear, respectively, 

9 dVt 

If = If int )p q ; (4) 

where fi^r) is the ion toroidal angular frequency, and r = o) is the radial 

coordinate that labels flux surfaces. Here, a and i \) 0 are the minor radius and the poloidal 
flux that label the last closed flux surface. The advective term has the momentum pinch 
coefficient /T,, and the diffusive term is proportional to the momentum diffusivity Xtp- 
Here, Ti int is the intrinsic momentum transport, also called the residual stress, which 
is the momentum flux generated even for zero flow and flow shear (i.e. = 0 and 

di'Ip/dr = 0). In steady state, the total toroidal angular momentum transport should 
be zero if there is no external source (i.e. n = 0 in equation Q), and the balance 
between non-zero pieces in n determines the radial profiles of the toroidal flow, fi^r). 

The momentum advection and diffusion have been theoretically investigated in 
previous work [131 13 (IS US], and theory has been compared with experimental 
observations in many tokamaks with strong ExB flow (Mach number ~ 1) driven by 
strong neutral beams on m, na ng. The momentum advection is not just due to 
the momentum carried out by the particle transport, but also due to inherent inward 
momentum transport (momentum pinch). The momentum pinch is driven by the drift 
due to the Coriolis force in the rotating frame [15]. The flow shear is diffused out by 
the same mechanism by which the ion thermal energy is transported. Consequently, 
the diffusion rates for the momentum and the ion thermal energy are similar, and the 
Prandtl number (the ratio of the momentum diffusivity to the ion heat diffusivity) is 
found to be 0.5-0.8 in many gyrokinetic simulations 0 E m eh. We have found that 
this previous analysis for and x<p corresponds only to the ExB flow, and it is not 
valid for the diamagnetic component of the neoclassical flow [221123]. 

The intrinsic momentum transport accounts for all other momentum transport 
mechanisms that are not the pinch or diffusion. This intrinsic momentum transport 
determines the intrinsic rotation. Many features of experimental profiles of rotation 
(e.g. the sign change of the rotation at mid-radius in some discharges) cannot be 
explained by only momentum pinch and diffusion, and so require intrinsic momentum 
redistribution. The intrinsic momentum transport must be driven by the mechanisms 
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that break the symmetry of the turbulent momentum flux in a non-rotating state (e.g. 
up-down asymmetric magnetic equilibrium [24, [25], the slow radial variation of plasma 
parameters [2.6, 27] and the slow poloidal variation of turbulence [29], neoclassical flows 
[231130113T1 [32] , hnite orbit widths [30], [31]). Reference [8] treats all these symmetry 
breaking effects in a self-consistent way. 

In this paper, we focus on the intrinsic turbulent momentum transport due to 
neoclassical phenomena. We find that the diamagnetic component of the neoclassical 
parallel particle flow and the ExB flow break the symmetry of the turbulence differently 
[23] . As a result, it is possible to obtain turbulent momentum transport for non-rotating 
states in which the diamagnetic flow and the ExB flow cancel each other. Moreover, the 
neoclassical parallel heat flow also results in significant toroidal momentum transport. 
Because the size of the diamagnetic flow (Mach number < 0.1) is smaller than the ion 
thermal velocity by (B / Bg)pi/L t, calculating the effect of the diamagnetic flow requires 
higher order corrections to the gyrokinetic equations [22], |23[ [33] - Here, B and Bg are the 
magnitude of total and poloidal magnetic field, respectively, and Lt is the characteristic 
length of the temperature gradient. The momentum flux due to diamagnetic effects may 
explain the collisional dependence of the observed intrinsic rotation PE21E2] and the 
rotation peaking during L-H transitions [23]. We investigate the effect of each piece of 
the diamagnetic corrections separately in section 4. 

The rest of this paper is organized as follows: In section 2, we introduce gyrokinetics 
and review the symmetry of the turbulence that results from the lowest order gyrokinetic 
analysis. In section 3, the neoclassical particle and heat flows are discussed. They 
are used to correct the lowest order gyrokinetic equation. In section 4, using higher 
order corrections, we evaluate the intrinsic turbulent momentum transport due to 
neoclassical effects for low and moderate collisionality. We show the momentum flux has 
contributions from (1) the neoclassical particle flow, (2) the neoclassical parallel heat 
flow, and (3) the neoclassical poloidal electric field and poloidal gradients of density and 
temperature. Finally, a discussion is given in section 5. 

2. Gyrokinetics and Symmetry 

Gyrokinetics is a kinetic description of electromagnetic (or electrostatic) turbulence 
that averages over the gyromotion while keeping hnite Larmor radius effects Em EE. 
The frequency of drift-wave microturbulence is much smaller than the gyrofrequency, 
permitting the separation of the gyromotion from the parallel streaming and the 
perpendicular drift motion. Particle motion is described by the position of the gyration 
center, and only two velocity space variables (for example, kinetic energy £ = v 2 / 2 
and magnetic moment p = v\/2B). Here, v± = ^Jv 2 — ujj and uy are the velocities 
perpendicular and parallel to the static magnetic held, B, respectively. Gyrokinetics 
assumes pi/l± ~ 1 and l±/l\\ ~ Pi/a = p* « 1. Here, l\\ and l± are the turbulence 
length scales parallel and perpendicular to the static magnetic held, respectively. For 
simplicity, we will only use gyrokinetic equations for electrostatic turbulence which is 
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the most relevant in low (3 tokamaks. 

The radial flux of ion toroidal angular momentum in equation ([2]) is dominantly 
due to the fluctuating radial ExB drift carrying toroidal angular momentum, and it can 
be evaluated using 

n ~ Iff 6 = (mi J d?vfl b (v ■ <pR)(v% ■ , (5) 

where • Vff is the radial drift due to the electrostatic turbulence and f/ b is the 
turbulent piece of the ion distribution function. The fluctuating radial ExB drift and 
distribution function are given by the gyrokinetic equations. The radial flux of the 
toroidal angular momentum can be also decomposed into a piece due to the fluctuation 
of the parallel velocity and a piece due to the fluctuation in the perpendicular velocity, 
Iff 6 = II|| + II_l. The parallel contribution to IT 6 in equation ([5j is 

n, = (m t I d?vff (vf ■ V*)) , (6) 


and it dominates over the perpendicular contribution II ^ ~ (Bg/B)Tl\\ when Be/B -C 1. 
In most tokamaks, the poloidal magnetic field Bg is much smaller than the total magnetic 
field B. 

The radial momentum transport can be expanded in the small parameter p*^ = 
(B/Bg)p+ that measures the radial drift orbit width with respect to the size of the 
tokamak, i.e. Iff 6 = IIi+n 2 + .... The lowest order radial flux (Eli) corresponds to the size 
of the gyro-Bohm diffusion of thermal velocity, hi! ~ D gBi {nimiV t iR/a) = plnpriiv/iR, 
where D gB ^ = p*PiVu is the gyro-Bohm diffusion coefficient. The next order radial flux 
(n 2 ~ (B/Bg)p± IIi) is given by the higher order gyrokinetics described in section 3.1 


The lowest order gyrokinetic equations described in this section can be used to evaluate 
the gyro-Bohm scale transport of particles and energy, but the lowest order radial flux 
of momentum (Iff) is different from zero only in the presence of sonic plasma flow, 
otherwise IIi = 0 as we will prove. 

The distribution function of species s is assumed to have a non-fluctuating 
background piece f bg and a fluctuating piece due to turbulence f tb (i.e. f s = f bg + f tb ). 
To lowest order in p*, the distribution function is approximately f s ~ /q® + f{ b s . 
The subscript denotes the order in p*. The plasma is sufficiently collisional that the 
background piece is /q® = /m,s, where f m,s is a stationary Maxwellian. The size 
of the fluctuating piece is much smaller than the background piece, f{ b s ~ P*/m,s- 
The electrostatic potential can also be split into a background and a fluctuating piece 
<f) = (j )/ 9 + (f) tb , where (f ) b 0 9 is the background potential, of order T e /e, (pf is the short 
wavelength turbulent potential whose size is 0(p*(T e /e)), T e is the electron temperature 
and e is the electron charge. The lowest order gyrokinetic equation for fl b s is 
ctb 


dfts 

dt 


k b 9( 


- x b) ■ v/“ + (t>,b + v„,- -V(y s > x b) ■ v/“ 


= yv{0“> X b ■ («|b + \ m ,,) ■ + <c(/,)>, 


Z,e 


B 


m x 


B 

BfZ 


ctb 


d£ 


(7) 
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where C is the collision operator and vm, s = (/z/f2 s )b x VB + (uy/h2 s )b x (b ■ Vb) is 
the sum of the VB and curvature drift. Here, m s is the mass and Z s e is the charge 
of the species, = Z s eB/m s is the gyrofrequency, and (...) is the average over the 
gyromotion. Note that (p b 0 9 is evaluated at the flux surface of interest ip = ipo in the 
second term on left hand side of this equation (Jr]) . We ignore the dependence of <p b Q 9 on 
ip to this order. 

The gyrokinetic equations for different species are coupled by imposing the 
quasineutrality condition, 



Z a e(<pf - 
T s 



0, 


( 8 ) 


where the second term in the parentheses on the left hand side is the polarization density 
due to the gyromotion. 

Because the turbulence parallel length scales are much larger than the perpendicular 
length scales (l±/l\\ ~ p*), it is convenient to use the field-aligned coordinates (ip, a, 
and 6) to solve equation ([7]) , with the poloidal angle 6 measuring the distance along B. 
The coordinates ip and a are perpendicular to the magnetic held. The angle a is 


a = ip — I(ip) / d0' 


R 2 B • VO' 




(9) 


The static magnetic held satishes B = Va x Vip. The slow variation in the parallel 
coordinate (9) and the fast variation in the perpendicular coordinates (ip and a) of the 
turbulence are well described by solutions of the form 

f*W,a,0-,£,ii,<r,t) = J2 kc.,V;£,H, ( 10 ) 

kipikot. 

= Y. ’ (U) 

kip ,/Cq: 


where and k a are the wavevectors in ip and a coordinates, respectively, and a is the 
sign of the parallel velocity (u|| = a\j2(£ — /jlB )). 

For a poloidally up-down symmetric tokamak, the lowest order gyrokinetic 
equations ([7]) and (|8|) satisfy the symmetry mm 


fi b s (6, v \\,k^) = -fi b s (-6, -V\\,-k^), 

4>i (9,k^) = -4>?(-9,-k^). 


( 12 ) 

(13) 


To see the ehect of this symmetry on the transport of ion toroidal angular 
momentum, the parallel contribution in (|6]) is written as an integral and a summation, 


n || = / d0dv\\n\\(9,v\i,k^), 


(14) 


where 


Try 


27T rriil 2n 

B ■ V9 V' 


Y / dn(ft b Yv\\ik ct <t> th jQ(k ± .v±/Sl i ). 


(15) 
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The asterisk on f* b indicates complex conjugation, Jo is a Bessel function of the first 
kind, and k± = \J k^V+ k^\S7a\ 2 . The parallel contribution IIy and the integrand 
7r|| are expanded in p^o (i.e. fl|| = np + np-H..). Accordingly, ri|| x = Yh ki f^o d6dv\\TT\\\ 
is the lowest order radial flux and II ||2 = d9d v \\TTp is the higher order flux in 


P*,e- 

From the symmetry properties in equations ( 12p3 ), the momentum transport is 
antisymmetric under the inversion of the poloidal angle, the parallel velocity, and the 
radial wave vector, 


7ij|i(0,t;||,/ty) = -7 T||i(-6>, -U||, —kp). (16) 

The perpendicular contribution II^ has the same symmetry as the parallel contribution 
IIj| has mm- As a result, the turbulent momentum transport vanishes (i.e. kb = 0) 
due to the symmetry of the gyrokinetic solutions unless there is a symmetry breaking 
mechanism. 


3. Symmetry breaking 


There are several effects that break the symmetry of the lowest order gyrokinetic 
equations described in the previous section. If the toroidal rotation is large (Mach 
number ~ 1), the symmetry of the gyrokinetic equation is broken even at lowest order 
(i.e. the symmetry in (16) is not valid because (p b 0 g is not 0(T/e), but 0{p~ l T/e)). 
However, we are interested in describing intrinsic rotation, which is typically small 
(Mach number ~ 0.1). The symmetry breaking mechanisms for the intrinsic toroidal 
rotation are described below. 

The symmetry in 6 is broken in up-down asymmetric tokamaks m- Typically the 
up-down asymmetry is significant only around the edge. Some dedicated experiments in 
the TCV tokamak m with significant up-down asymmetry in the core found rotation 
of about 3 % of the ion thermal speed. A theoretical study in [25j estimated rotation 
due to the asymmetry to be of the order of 5% of the ion thermal speed. 

The slow variation of plasma parameters in the radial direction breaks the symmetry 
in kp HU \ZE\- The effect on the turbulence of the change of temperature and 
density gradients across the radial dimension of plasma eddies has been investigated 
in gyrokinetic global codes. For a sufficiently small radial correlation length of the 
turbulence, the effect of the slow radial variation of the gradients is too small to break 
the symmetry significantly. 

The higher order correction to the gyrokinetic equation due to the slow poloidal 
variation of the turbulence can break the symmetry in the poloidal angle 6 [29]. This 
effect will be small if the parallel correlation length of the eddies is small compared to 
the connection length qR [8] , with q the safety factor. 

In this paper, we focus on the effect of the neoclassical flow that introduces a 
preferential direction in the system [22| [23], [321133] . For turbulence in which the radial 
and parallel correlation lengths of the turbulence are sufficiently small, these diamagnetic 






Figure 1. Sketch of the origin of the diamagnetic particle flow. The poloidal cross section 
of a flux surface is represented by the black solid circle. The particles with positive parallel 
velocity drift radially inward (blue circle), while the particles with the negative parallel velocity 
drift radially outward (red circle). Because the plasma at inner radii is denser and hotter than 
the plasma at outer radii, the total parallel flux at the outer mid plane is positive and the ion 
diamagnetic flow is ^ ~ Aif>(d Inpi/dtp)vti ~ p+{B / Bg)vu- 


flows dominate over the momentum driven by the slow radial and poloidal variations of 
the turbulence [8]. 

The toroidal flow in a tokamak is composed of two different types of flow: the ExB 
flow and the neoclassical or diamagnetic flow. The diamagnetic flow due to the radial 
pressure gradient always exists in a tokamak regardless of the size and the sign of the 
radial electric field. This implies that there is an inherent preferential direction that 
breaks the symmetry. We call these flows diamagnetic flows to emphasize that they are 
the result of finite orbit width effects, as shown in figure 1. The higher flux for the ions 
with positive parallel velocity is due to the radial magnetic drift to inner radii where 
the plasma has a higher temperature and density than at outer radii. Consequently, the 
amount of diamagnetic flow is determined by the width of the deviation of the poloidal 
orbit from the flux surface (~ (B / B e )pi ) and the pressure gradient (i.e. ~ P*,eVti, 

where is the toroidal angular frequency for diamagnetic flow). To take into account 
the diamagnetic flow in gyrokinetics, one needs to correct the lowest order gyrokinetic 
equation in ([?]) with the higher order terms in p*^ as will be explained in section 3.1 


3.1. Higher order gyrokinetic equations 


The symmetry in (12) and (13) is broken by the effect of higher order terms in p*^. One 


thus needs the higher order gyrokinetic equations to evaluate the momentum transport 
(i.e. hi! = 0 and n 2 ^ 0). 

In the gyrokinetic equations in this section, we only keep the corrections related to 
the diamagnetic flow and neglect all other corrections in p*^ for simplicity. As shown 
in p|, they are the dominant terms when the perpendicular characteristic length of the 
turbulence l± is small compared to the characteristic poloidal gyroradius (B / B e )pi . The 
higher order ion gyrokinetic equation with both the diamagnetic flow and the ExB flow 
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in the lab frame is 

dftb c 


dt 


c 

B 


dp) 


W’o) + (ip - V>o) 


dil ; 2 


(V>o 


v^xb -V/; 


ctb 


+ (n||b + v M - —V((p b l 9 + ( <p tb )) x b 


•v/ f 


tb 


(f||b + vj,)-V^ 
TTlj ac 

+ ^(n,|b + v M )- V(0 t6 > 9/< 
rrii 


dS 


= -V(0“) x b • V/, 
9 + <<?(/*)>, 


t>9 


(17) 


where the distribution functions and potentials include the higher order corrections; i.e. 


ftb _ ftb , £tb fbg _ cog . rug. rot 

Ji ~ Jl,i "T" 12,j) Ji — J0,i ' Jl,i ' J 1, 

that in the second term on the left hand side of the equation, we have Taylor expanded 

d(p b 0 9 /dip to keep the dependence on ip. Here, /q® = /m,i is the lowest order background 
piece that is a Maxwellian distribution function as in (7), and /^’ E and f\ 9 ) d are the 
small deviations from the Maxwellian due to the ExB flow and the diamagnetic flow, 
respectively. For intrinsic rotation, the ExB flow is of the same order as the diamagnetic 
flow (i.e. f\ 9 p E ~ f\ 9 / 1 ~ P*,6>/o~D- For the ExB flow, the correction to the background 
distribution function simply results in a shift of the background Maxwellian by the ExB 
flow (i.e. /J? + f\ 9 f = f M ,i( v - %,eR<p)), making f\ 9 f ~ (m i U||/T i )(/n Vi£ -/S)/ Jtfii 
where = — c(d(pQ 9 (ip) / dip). The small deviation due to the diamagnetic flow, /^’ d , 

is the neoclassical distribution function Esa eu, which will be discussed in the next 
section. The higher order neoclassical potential <p b9 [ip , 9 ), which is determined by (23), 


cbg,E 


cbg,d 


= $ + 


k b 9 


and (p tb = (pf + <pf . Note 


is also considered in the above gyrokinetic equation (|17|). 
The quasineutrality condition in 


must be corrected to include the higher order 
pieces and <pf, but the neoclassical corrections to quasineutrality are negligible for 
Bg/B <C Pi/l_ l <C 1 (see [8] for the details). We consider this limit, thus can use (|8|). 

From the size of the correction to the background piece, we obtain the size of the 
corresponding corrections to the fluctuating pieces, ~ P*,eft b -i and (pf ~ p*,o<p tb - For 
the transformations uy — > — uy, 9 —y —9 and —> — k^, all additional higher order terms 
in the gyrokinetic equation have the parity opposite to the lowest order terms in ([7]). 
The parity of the higher order corrections to the distribution function and potential is 
obtained by linearizing the effect of the corrections small in p*.o (compare the following 


with (12) and (13)), 

/?i(M||,*ty) = -I’ll, -M (18) 

$(9,k+) = $(-9,-^). (19) 

The different symmetry of ff i and <pf results in non-vanishing higher order momentum 
transport, 

n\\2{0,V\\,ki,) = 7Ty 2 ( — 9, -uy, —/ty), 
where 7Ty 2 oc Re[ik a (f*i<f>¥ + f*<pf )^||]- 


( 20 ) 
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3.2. Neoclassical correction to Maxwellian equilibria 


The correction f\ 9 f to the Maxwellian distribution function in (17) is derived from the 
drift kinetic equation for species s !Mi m, 


U||b • + v MjS ■ V/q 9 = C 8 {H 1>8 ), 


( 21 ) 


where H ] s = f\ 9 f + (Z s e0//T s )/ 0 9 . Here, the lowest order background potential 0 O 9 
which results in the ExB flow is not included in (21), but the higher order background 
potential qb\ 9 due to the neoclassical correction is included. The lowest order potential 0 q 9 
is considered separately in f\ g J. R of (17). Using v M , s • V/q 9 = Iv yb • V(u||/fU)(<9/o 9 /dip), 
the drift kinetic equation is simplified to 

»||b • V - /") = C,(H U ), (22) 

where /j'( — — (Iv\\/£l s ) [ dhin,/dTp + (m,-v 2 /2T s — 3/2)(c*In T,/chl>)\ Depending on 
collisionality, different solutions for Hi s are obtained. By employing the quasineutrality 
condition, 

Z 8 


Y,z.jd\,U„- z ^f b A=a, 


(23) 


we can also obtain (p\ 9 . 

To study the symmetry breaking mechanisms due to the higher order background 
distribution function, the solution of the ion drift kinetic equation is decomposed into 


fbg,d _ 1-|| poy^ii og 

Jl,i ~ Jl,i "r Jl,i "r" Jl,i 


9,V\\ 


?ll 


pbg, other 




(24) 


where /j^ 1 is the piece resulting in the diamagnetic parallel particle flow, /^’ 9|1 is 
the piece giving the diamagnetic parallel heat flow, and f b9 .’ other accounts for all other 
contributions. 


1 he piece j 1 i is 
fb9,V\ | _ m N\\V\\ 4 fh g 

J 1,1 rj -1 J 0 


(25) 


where V\\ id {dpi/dip, dT t /dip, u*...) is the diamagnetic parallel particle flow that depends 
on the pressure gradient, temperature gradient and collisionality z/*. Here, z/* = 
(unqR/e 3 ^ 2 Vti) is the ratio of effective collision frequency for pitch angle scattering to the 
bounce frequency, u vi is the ion-ion collision frequency, and e = r/R is the inverse of the 
aspect ratio. To lowest order, for a plasma with only one ion species, the diamagnetic 
parallel particle flow is given by 


V\A - fw) - K Mb, 




(26) 


where = —(c/Zieni)(dpi/&ip) and K{ip) oc dR/dip. For example, the diamagnetic 
parallel particle flow in the banana regime zz*Cl for concentric circular flux surfaces 
and large aspect ratio is 


V ld = 


it .;• 




1 d In pi 
h dip 


1-17 f c h 


9 In Ti 
dip 


(27) 
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where h — B/Bq, Bo and hho are the magnetic held magnitude and ion gyrofrequency 
at the magnetic axis, respectively, and f c is the effective fraction of passing particles 

m- 

Combining the parallel diamagnetic flow with the parallel flow due to the lowest 
order radial electric field and the perpendicular flow, the total equilibrium flow in a 
tokamak is obtained [38], 

Vi = (Vf| >d + V|| jjE )b + V_L 

= + %,eW)Rv + (28) 

where V^e = — (cl/B)(d(po 9 /d't/j) is the parallel flow due to the lowest order radial 
electric field. Here, Vi = — (f^, p + O^#) (B x V'i/j/B 2 ) is the perpendicular flow due to 
the ExB and the pressure driven diamagnetic drift, and we use b x V'tp = I b — BR(p. 
Notice that the poloidal flow is determined only by the term proportional to K (?/), while 
the toroidal flow is given by £ltp,p + &tp,E and the term proportional to Thus the 

toroidal particle flow is 


V v = Vi-<p = %R 

= (f\,,p(V0 + %,e{^) + 0))R (29) 

where /{n^^R 2 ^, 6)). We can divide the toroidal angular 

frequency of the temperature gradient driven flow into two pieces, i t('0, $) = 
^,tC0) + AO ^(^,9). The piece contributing to the angular momentum is 


(!J„, T ,R 2 )*(V>) _ KW)I(4 ,) CdB( B ■ V#)- 1 




(R 2 )4^) rii( i/j) J^ddR 2 (B-Vd)- 

and the piece not contributing to the toroidal angular momentum is 


\-i 


A = 


UiW \R 2 ^,0) J o 2n dOR 2 (B-VO)~ 1 


Jo d0(B-vey 

2 TT 


where (AQ^^R 2 )^, = 0. Then, the particle flow is given by 

V v = + %,eW)R + AHp, T (^, Q)R 


(30) 


(31) 


(32) 


where = f^ )P (^) + 

The diamagnetic parallel heat flow is due to the terms proportional to v\\v 2 /2 in 
the distribution function, The parallel heat flow breaks the symmetry even in the 

absence of rotation because it corresponds to a distribution function that is odd in dm, 


bg , qi | _ 2miV\\q\\ f rriiV 


r'-'b/M || ~ 

Jl.i r 


5 PiR 


2 T 


^) fb 


o > 


(33) 


where pt is the ion pressure and the heat flow qn (dTi/dip,^...) is also a function of 
the temperature gradient and collisionality. We can see that (25) and (33) satisfy 
J dv \\ v l\fl 9 i V " = V ld, I dv\\v\\f^- 911 = 0, J dv ]] v\\(m i v 2 /2T i - 5/2) = 0 and 

f dv ll v ll (m i v 2 /2T i - 5/2) /fj’ 9 " = q l{ . 
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3.3. Symmetry breaking due to diamagnetic effects 

The two types of toroidal rotation and ^L Vi e) have different origins and 

characteristics. The rotation H v .e is proportional to the radial electric field because 
of radial force balance, and as a result the radial electric field changes in the momentum 
confinement time scale The electric field changes the energy and orbits of each 

particle. Conversely, the pressure and the temperature gradients in iip.d are determined 
by the turbulent anomalous transport of particles and energy, changing on the energy 
confinement time scale. The pressure and the temperature cannot change the energy 
and orbits of individual particles. These different characteristics are reflected in the 
gyrokinetic equation ( [l7| ) by the different terms for and There are additional 

terms on the left hand side of (17) related to (compare to ([ 7 ])). The term 

—(c/B)(f> — fj 0 )(d 2 cj) b 0 9 / df> 2 ) (\7 x b) ■ Vff b is the effect of the background ExB drift. 
The other term, — [Zie/mfffvM • V(p b 0 9 )(df bb /dS), is the acceleration of the particles 
in the radial electric field when the particles change their radial positions due to the 
curvature and VB drift. The diamagnetic flow does not have terms like these on the 
left hand side because it cannot change the orbits and energy of the particles. 

Due to their different characteristics, the effects of the two flows and on 
the turbulent momentum transport are different. Then, in the low flow regime, it is 
convenient to write the ion toroidal momentum transport as 


11 u int TiiTriiiB. P<p,E^tp,E') nirrii{R P^^AXt^pdfip 


_ , n m l f? 2 \ Iv . SO. lf E 

| Xp,d "T X<p,E 


dr 


- npmi ( R XpM 


dAQ 


<P,T 


dr 


(34) 


where IT( nt is the redefined intrinsic momentum flux for = 0, = 0, A= 0, 


d^lp,d/dr = 0, dCl^^/dr = 0, and dAfl^^/dr = 0. Equation (34) is valid for f 


and sufficiently small to be able to linearize II around = A = 0. 

By using ( R 2 Xl ip ) 1 p = (R 2 )qp(Q<p,d + D^e), the original model in (jfj) is recovered by 

II H int nimi(R ) qpPp^E T ^p,e') 

/ p2\ f dQj^d 9 Q<pe\ 

- niiTiifR }qpXp,E ( H-yy— ] ■ (35) 


dr 


dr 


Then, equations (34) and (35) result in the relation between Yl int and II' ni 

TT — TT' 

LL int 


int 


Tl%1Tli (.R ^ip/\pTliTfli (-R Pip 


- nim^R 2 )^ Ax, 


dn 


V,d 


dr 


— nririi 


R 2 Xv, 


dAn 


‘ P,T 


dr 


(36) 


where A P^ = P v4 ~ Pp,E and Ax v = X^d ~ X<p,e ■ Here, lT mt includes the effect of the 
parallel heat flux, 1 , of f\ 9 f othe ' in (24), and of the neoclassical potential (f>\ 9 . 


The intrinsic momentum transport in (|36|) is then 

+ IT 


TT — TT A ' F 7 -4- TT Ai<:¥ 

11 — LL int + i b 


int 


l e 

int 


+ n 9 " 


int 


1 T\^i 1 Toother 
' L1 int "r iA * 


int 


(37) 


where the contribution of the difference in the momentum pinches is 


n,tf* = 


(38) 
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the contribution of the different diffusivities is 


n fnt* = -ni7m(R 2 )^ 


on 


ip,d 


dr 


(39) 


and the contribution of the different poloidal dependence of the temperature driven 
diamagnetic flow is 

dAfl Vt T 


n 


int 


Here, and n Z are the intrinsic momentum transport due to the symmetry 

breaking caused by the parallel heat flow and the neoclassical potential, satisfying 

ll) and ^li 9 ) = respectively. Finally, Iif^ t er 

includes the contribution due to the piece f^ 9 .’ other breaking the symmetry of the 


= -riim^R 2 P^bAVL^t)^ ~ n,jn t (r 2 x v ,6 


dr 


(40) 




±bg 


turbulence, because it satisfies /, er (6,v\\) = —fii° er (—6,—v\\). 

Thus, the intrinsic momentum transport for zero flow and zero flow shear is a 
function of f l^d, dQ^j/dr, AO Vi t, Q\\ and p b9 . In general, the diamagnetic flow depends 
on the pressure gradient, the temperature gradient, and collisionality, so the intrinsic 
momentum transport is 


rr _rr l d Pi dT i 92 Pi ^ T i du * TV 

int ~ int 1 di) ’ di[> ’ I/ *’ d1P 2 ’ dip 2 ’ dip ’ int 


(41) 


In the next section, we will show the contribution of each term to the momentum flux 
AP 

by evaluating 11, v 


int. 


LL int i 


TI° nt , n^ t , and with gyrokinetic simulations. 


4. Evaluation of intrinsic momentum transport 


We used the gyrokinetic code GS2 [39] to evaluate the intrinsic momentum flux. For 
numerical reasons, the gyrokinetic equation is solved in the frame rotating with f hg , 


thereby avoiding the energy derivative terms due to p b 0 9 on the left hand side of (17) 
The kinetic equation in the frame rotating with 0 ) 

MR) 


on 


dt’ 


+ w - ■ v/f (s) 


dip 


v'\\b + V M + v( 


C,£1 V: e 
bg(R) 


g V(^ + (0“» X b 


v/, 


t.b(R) 


—vUWPI 9 

mi 11 


= — V (p tb ) x b ■ V/, 

+ — (w||b + v' M + v 7 ^ 
rrii " 


+ -V(0 tb ) x b • Vip mJVl1 dn j’ E fWV 
B w ' T t B dip Jl 


V(0“) 


Qfbg(R) 


+ <■ C(fi )), 


o MR) 

b)—- 

; d£’ 


(42) 


' \~r / Qgi 

where d/dt' = d/dt + )RpP • V is the Doppler shifted time derivative, E = 

(2uy Q^^/Q^b x [(V/? x <£) x b] is the Coriolis drift and f b9( ' B ' 1 = /o,^( r; f|) + fip\ v \\) 
is the perturbed Maxwellian distribution function due to the diamagnetic effects in the 
rotating frame. Here, we apply a transformation from the lab frame velocity v to the 
rotating frame velocity v 7 using the relation v 7 = v — The parallel velocity 
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becomes v'» = v\\ — I£1 Vj e/B and the kinetic energy becomes S' = £ — I£L Vj eV\\/B 
[31]. The superscript (R) indicates the value in the rotating frame, whereas no 
superscript denotes the value in the lab frame. The conversion of the gyrokinetic 
equation between the laboratory frame and the rotating frame is given in appendix 
A of [31] for 

To investigate the intrinsic momentum transport in a non-rotating state, we set the 
toroidal flow to be zero, £1^ = + £ltp,E — 0- The intrinsic momentum transport is 

the flux in the lab frame, which is 

Hint = H tb = n tb ( R) + min v , E {TR% = n tb(R) - m^RTR 2 )^ (43) 

where T is the particle flux. 

For the gyrokinetic simulations in this paper, we used two species (deuterons 
and electrons), and the plasma parameters of the Cyclone Base Case that model 
experimental conditions in the plasma core : Rq/Lt = 6.9, Ro/L n = 2.2, q = 1.4, 
r/a = 0.54, R 0 /a = 3.0 and s = 0.8. Here, Lt = — Ti/(dTi/dr ) = — T e /(dT e /dr ) and 
L n = —nj ( dn/dr ) are the characteristic lengths of temperature and density, respectively, 
R 0 ~ {R)t/j, q is the safety factor, and s is the magnetic shear. Due to the new higher 
order terms in (42), we need to specify values for p*, d(R 0 /L n )/dr and d(R 0 /L T )/dr. 
We choose p* = O.Olf, d(Ro/L n )/dr = 0 and d(Ro/ Lt) /dr = 0. 

The GS2 simulations use 32 grid points in the parallel coordinate 9 , 12 grid points 
in kinetic energy, and 20 grid points in pitch angle. The box size of the simulation in 
both radial (-0) and binormal ( a ) directions is approximately 125 p*. We use 128 and 
22 Fourier modes in the radial and binormal directions, respectively. 

Neoclassical effects depend strongly on collisionality z/*. We investigate the low 
collisionality intrinsic momentum transport due to neoclassical particle and heat flows 
in section 4.1. We compare these results to results for moderate collisionality in 
section 4.2. In the moderate collisionality case, the contributions of f^ 9 .’ other anc ) 
to the momentum transport are non-negligible, whereas they are not important for low 
collisionality. 


f In many experiments, 0.001 < p* < 0.01 
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(a) 



9 [Radian] 


(b) 



Figure 2. (a) The normalized diamagnetic parallel particle flow (red) V^d/vu and the 
parallel heat flow q\\/vuPi (green) in terms of poloidal angle, (b) The normalized radial 
gradient of diamagnetic parallel particle flow (red) (dV\\,d/9 r )/( v ti/ a ) and the parallel heat 
flow (dq\\/dr)/(vtiPi/a) (green) in terms of poloidal angle. All profiles are obtained with the 
neoclassical code NEO f40] [4T] for the Cyclone Base Case with p* = 0.01, d(Ro/L n )/dr = 0, 
d(Ro/LT)/dr = 0 and low collisionality = 0.043. 


4-1. Low collisionality 

Table 1 shows GS2 results for the contributions of each neoclassical correction to the 
momentum flux for low collisionality z/* <C 1. For all the contributions, the normalized 
radial flux of toroidal angular momentum divided by the normalized turbulent heat flux 
is of order the poloidal rhostar {B/B e )p i , ~ 0.1, as expected. The neoclassical parallel 
particle and heat flows and their radial gradients in figure 2 are used to obtain the 
simulation results in Table 1. In the low collisionality case, the total momentum flux is 
almost the summation of the contributions from the particles and heat flows, which are 
comparable with opposite signs (i.e. II int ~ IlJ t + II^ t from Table l-(a), (b) and (c)). 
Other contributions to the momentum flux in Table l-(d) and (e) are clearly smaller. 

4-1.1. Neoclassical parallel particle flow contribution. The momentum transport due 
to the neoclassical parallel particle flow can be divided into three pieces: (1) flux 
due to the different pinches for the diamagnetic flow and the ExB flow (Llf^ = 
—7iimi(R 2 )pAP ip fl lfj d), (2) flux due to the different momentum diffusivities for the 
diamagnetic flow and the ExB flow (II= —n i m i (i? 2 )^A / \' ¥ ,(cIf2 v , )(i /c)r)), and (3) flux 
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Table 1. The intrinsic radial flux of toroidal angular momentum II (we assume a non-rotating 
state = flip,d + Q<p,E = 0), divided by the ion heat flux Qi given by GS2 for the Cyclone Base 
Case with p* = 0.01, d(R 0 /L n )/dr = 0, d(R 0 /L T )/dr = 0 and low collisionality v+ = 0.043 for 
(a) full neoclassical function and potential given by NEO, (b) only the particle flow piece , 

(c) only the parallel heat flow piece // 9 ’ 911 given by NEO, (d) only the neoclassical potential 
and (e) only the remaining piece f ^ other . Here, a is the standard deviation of the 

normalized momentum flux from the time averaged value due to turbulent fluctuations. 


Size 


GS2 results 


IEnt Vti 
Qi Ro 



(a) 

II int 

due to f/ 9 - d and (j)\ 9 

nbg,d _ rbg,V \|| 

Jl,i ~ Jl,i 

, f b g,q || 

Jl,i 

■ rbg, other 

' Jl,i 

-0.022 

0.018 

(b) 

ni 

due to f/ 9 - V " 

Vj|,d and dV|| )d 

r /dr in figure 2 (red) 

0.025 

0.020 

(c) 

n 9 " 

due to f\ 9 ; n 

q\\ and dq\\/dr 

in figure 2 (green) 

-0.046 

0.015 

(d) 

LL int 

due to cf>\ 9 

e-\/((^ 9 ) 2 )v-/T =* 0.0005 

0.002 

0.019 

(e) 

n? n h p r due to f%° ther 

rbg,d rbg,V\ jj 

Jl,i J l,i 

pbg,q || 
J l,i 

and = 0 

0.001 

0.016 


due to the different poloidal dependence of the pressure gradient driven flow and the 
temperature gradient driven flow (Id e int ). The momentum transport due to the particle 
flow in Table l-(b) is approximately the summation of the three pieces in Table 
2-(a), (b) and (c). 

The pressure driven piece of the rotation is VL^^^Rq/vu = —(c/Zie7ii)(dpi/di/j)(Ro/v t i) 
(l/V2)(q/(r/R 0 ))p*(a/R 0 ) 2 (R 0 /L n +R 0 /L t ) = 0.167, as shown in the green curve of 
figure 3-(a). By solving the drift-kinetic equation in (21) for this case using the neoclas¬ 
sical code NEO [40, H], we can obtain the piece of the particle flow that is proportional 
to the temperature gradient, (K (-0)/ n^))B in (26), as shown in the blue curve of 
figure 3-(a). As shown in figure 3-(a), this piece tends to cancel the other piece of the 
particle flow fl<^ p (^)I% in a low collisionality regime, giving Q^^Rq/vu = 0.087 and the 
strong poloidal variation Ain figure 3-(b). The red curves in figure 3-(a) and (b) 
correspond to the red curve in figure 2-(a), representing the particle flow V» d in which 
both contributions from and A fltp,T are included. 

The different pinch coefficients were investigated in [23] for the case with Rq/Lt = 

9.0, Ro/L n = 9.0, q = 2.5. For these parameters, there is a 22% difference between the 
pinch coefficients: P v ,d/x<f,E — 3.5/1% and P^e/Xv^ — 2.9/1%. This difference in the 
pinch coefficients may explain the rotation peaking observed in H-mode [23]. For the 
Cyclone Base Case, the difference between the pinch coefficients is also non-negligible, 


A P 


^ _ *<p,d 


P 


<p,E 


Ro/Lq 


X<p,E 


Xq>,E 


Qtp,dRo/^ti P^E 



0.3 

w 


(44) 


where we used the Prandtl number Pr E = x^,e/Xi — 0.54, which is obtained for the 
case (<9fl^£/<9r)(I%a/u f i) = 0.152. The time averaged momentum flux due to the 
different pinches IIis obtained from a GS2 simulation with Q^^Ro/vu = 0.087 
and = —0.087 while eliminating contributions of other neoclassical effects 
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(i.e. f^ d = and A= d£l Vt d/dr = dXl VjE /dr = (p^ 9 = 0). The small size of 

the neoclassical particle flow for the Cyclone Base Case results in the relatively small 
momentum flux due to the different pinches, compared to the case in [23j. 

The origin of the different pinches is the acceleration term that depends on the type 
of rotation, (Z i e/mic)Q tp>E (vM ■ Vip)(df* b /d£) in (17). In the presence of ExB flow, the 
particle is accelerated within the radial electric field, while the diamagnetic flow does 
not give acceleration. The acceleration term breaks the symmetry of the turbulence 
because it is odd in the poloidal coordinate (#). For example, for a circular tokamak, 
the term is proportional to sin#. 

The different momentum diffusivities result in the piece of the intrinsic momentum 
transport II , shown in Table 2-(b). The size of the momentum flux is proportional to 
the size of the velocity shear, which is determined by the second derivative of the pressure 
and temperature in radius. We set the size of the velocity shear dr)(R 0 a/vu) = 

—0.152 by setting d(R 0 /L n )/dr = 0 and d(R 0 /L T )/dr = 0 as done in [22]. The 
different momentum diffusivities are evaluated by comparing the Prandtl number for 
the diamagnetic flow and ExB flow, 

Xtp,E (Ro/L t) i 


APr = 


ttWv 

Writ v ti 


Xi 


(dfl V}d /dr)(R 0 a/v ti ) R 0 \ Qi 2 R 0 


~ 0.33, 


(45) 


where 


is obtained from a GS2 simulation with ( dXl lfi! d/dr)(R 0 a/vu ) = —0.152 
and ( dXl VjE /dr)(R 0 a/v ti ) = 0.152 while eliminating contributions of other neoclassical 
effects (i.e. f\ g f = and AXl VyT = = Q^e = (p\ 9 = 0). Because the Prandtl 

number for the ExB flow is about 0.54, this difference is significant, about 61%. 

The origin of the different diffusivities is the other rotation type dependent term, 
(1 /B)(ip — V’o ){d^,E/dip)Vip x b • V/* b , in equation (17). Generally, the radial shear 
of the rotation reduces the radial correlation length of the turbulent eddies because an 
eddy radially aligned at an initial time becomes misaligned in time due to the radial 
gradient of the toroidal flows 016]. The rotation type dependent term contributes to 
the radial misalignment differently for different types of rotation, because the particle 
orbits are modified not by the pressure gradient but by the radial electric field. We 
found that the difference of the Prandtl number becomes negligible (A Pr ~ 0.003), if 


tb 


the ExB shear term on the left hand side, (1 /B){\p — / ip 0 )(dXl iPtE /d'il’)'V^ x b • V/, 
is ignored but the velocity shear terms on the right hand side, (c/ B) (V {<p tb ) x b • 
'V'ip)(m i Iv^/T i B)(dXt iPtE / dip)fi 9 ^ for ExB flow and (c/B)(W {(p tb ) x b ■ V/^’ d ) for the 
diamagnetic flow, are kept in GS2 simulations. 

The difference in the momentum diffusivities depends significantly on the turbulent 
characteristics. For the plasma parameters in [23], the effect is negligible, which results 
in the small difference of the momentum diffusivities (i.e. A Pr/Pr < 0.01). 

The time averaged value of momentum flux due to a poloidal dependence H 9 int 
in Table 2-(c) is obtained from a GS2 simulation with the non-zero Aof the 
blue graph in figure 3-(b) while eliminating contributions of other neoclassical effects 
(i.e. f\ 9 / = /^ V " and = dXl^/dr = dXl^ E /dr = (p\ 9 = 0). This 
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Table 2. The normalized intrinsic momentum flux caused by the diamagnetic particle flow 
II^ t in GS2 for the Cyclone Base Case with p* = 0.01, d(Ro/L n )/dr = 0, d(Ro/LT)/dr = 0 
and low collsionality = 0.043. Here, has three contributions due to : (a) the different 
momentum pinches n^y'’, (b) the different momentum diffusivities and (c) the different 

poloidal dependence of the temperature and pressure gradient driven flows n® nt . 


Characteristics of non-Maxwellian f^ d GS2 results 



,dR() 

RqCL 

A Q cp ,t(4’, 0) 

IIi„t Vti 



Vti 

dr vti 

Qi Ro 

(a) ntf 

0.087 

0.0 

0.0 

-0.005 

0.019 

(b) nfj- 

0.0 

-0.152 

0.0 

0.043 

0.013 

(c) nL 

0.0 

0.0 

Fiq. 3-(b) blue graph 

-0.012 

0.015 


intrinsic momentum transport Iq nt occurs because the temperature driven diamagnetic 
flow 0) has a poloidal dependence different from that of the ExB flow f I v> e- 

According to the definition of Af^y^r, 6) in equation (31), there is no net toroidal 
angular momentum due to AQ^y^r, 0) on a flux surface. As shown in the blue graph 
of figure 3-(b), A VL Vi t > 0 at the outer-midplane (|d| < 7r), and AJI^t < 0 at the 
inner-midplane (|0| > n). The rotation AO^ can result in a finite flux of toroidal 
angular momentum because the ExB radial drift due to the turbulence is nonuniform 
poloidally. Usually, the ExB radial drift is larger at the outer-midplane than at the 
inner-midplane. Thus, there is more negative momentum pinch of positive momentum 
at the outer-midplane (—< 0 at 6 = 0) than positive momentum pinch of 
the negative momentum at the inner-midplane > 0 at 6 = n). On the 

other hand, there is more positive momentum diffusion due to negative radial gradient 
of the positive momentum at the outer-midplane (— x<p(dAQ V: T/dr) > 0 at 6 = 0) 
than negative momentum diffusion due to the positive radial gradient of the negative 
momentum at the inner-midplane (— XipidAQ^^/dr) < 0 at 6 = tt). Thus, the sign and 
size of n e int is determined by the two competing terms on the right hand side of (40). 











Vii B/l [Vti/Ro] 'cr 1 ' V|i [vti] 
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9 [Radian] 



9 [Radian] 


Figure 3. The normalized diamagnetic particle flow in terms of poloidal angle calculated by 
NEO for the Cyclone Base Case with p* = 0.01 and low collisionality z'* = 0.043: (a) The parallel 
particle flow (red) V\\^ = (Q V , P + the pressure-driven flow (green) and the 

temperature-driven flow (blue) = ( K/n)B , (b) The parallel particle flow divided by 

I /B (red) V\\^(B/I) = + An ¥ , i 7'(0), the poloidally uniform angular frequency piece (green) 

f2 , Pi d = + and the poloidally non-uniform angular frequency piece (blue) AQ^t^)). 

Here, Q v , p Ro/v t i = 0.167 and ( K/n)B 0 /v t i — —0.080 result in Cl Vt d.Ro/vu = 0.087 
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4-1.2. Neoclassical parallel heat flow contribution. The neoclassical parallel heat flow 
results in momentum transport II^[ t because it breaks the symmetry of the turbulence 
even in the absence of parallel momentum (i.e. J = 0). The piece for the 

parallel heat flow /^’ 9|1 is obtained by solving (21) with NEO. The normalized size of 


the heat flow q\\/v ti Pi is comparable to the normalized size of the particle flow V^/va, 
as shown in figure 2. The size of the heat flow is proportional to the temperature 
gradient and depends on collisionality, as the temperature gradient driven diamagnetic 
particle flow does. We found that the momentum transport due to the heat flow 

is significant, as shown in Table 1. To calculate the time averaged value of momentum 
flux due to the parallel heat flow IT^ t , we use a GS2 simulation with non-zero q\\ while 
eliminating contributions of other neoclassical effects (i.e. f\ 9 ? d = and <j) b fl = 0). 

For the cyclone case with low collisionality, IT^ (/ is negative and it dominates over fl] 


bi 

inti 


giving the negative total momentum flux II int . 

Figure [4] explains the physical origin of the momentum transport due to the radial 
variation of the parallel heat flow. Because the turbulent radial E x B flow Svexb has 
poloidal variation, the radial change of the parallel heat flow results in a temperature 
perturbation with poloidal dependence as shown in figure (4[(b). The temperature 
perturbation can be approximately described by 
Tnl 9 


ST, 


A 9q W 

SVExBTnl-?T- 
or 


(46) 


R„qde 

where r n i is the turbulent nonlinear correlation time. The pressure gradient in the 
parallel coordinate between the higher temperature location and the lower temperature 
location gives the perturbed parallel particle flow 

£ n i d6T i (A*7\ 

mM '" ~ T ^~W' (47) 

where the momentum equation in the parallel direction has been used. The product 
of perturbed parallel particle flow induced by the heat flow gradient and the turbulent 
radial E x B flow SvexB does not necessarily vanish in the poloidal angle integration, 
and it results in non-zero momentum flux IT^ t ~ ({mi5v\\8vExB)iii)T, as shown in figure 
[ 4 ]- (c). This mechanism gives IT^ t < 0 for dq\\/dr < 0, if huy at large poloidal angle is 
sufficiently small, which is consistent with the results in Table 3. 

Table 3 shows the simulation results of intrinsic momentum flux due to the parallel 
heat flows that are poloiclally uniform to remove the poloidal variation effect of the heat 
flows for simplicity. Significant momentum flux occurs due to the radial gradient of the 
heat flow, and the sign of the flux in Table 3-(c) is consistent with the mechanism in 
figure [4j explaining the result in Table l-(c). I 11 the absence of the radial variation and 
poloidal variation of the heat flow, the momentum flux is small, as shown in Table 3-(a). 


4-2. Moderate collisionality 


Table |4] shows GS2 results for the contribution of each neoclassical correction to the 
momentum flux for moderate collisionality z/* ~ 1. In the moderate collisionality case, 
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Figure 4. A sketch of a the mechanism that leads to momentum transport when there is a 
radial gradient of neoclassical parallel heat flow, li^[ t . (a) Assume that the neoclassical heat flow 
at inner radii is larger than the heat flow at outer radii. A positive heat flow without parallel 
particle flow means that hot particles move in the positive parallel direction and cold particles 
move in the negative parallel direction with the equal speeds, (b) The E x B flow Svexb, which 
is larger at the outer-midplane around 9 = 0 (ballooning structure), gives a poloidally varying 
perturbed heat flow Sq^. (c) The variation of the temperature in the parallel direction <57) ((9) 
due to the polloidally varying 8q\\ results in the poloidally varying parallel particle flow 5v\\(9), 
giving non-zero momentum flux as <5v||(0) is transported by the radial E x B drift Svexb■ 
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Table 3. Normalized intrinsic momentum flux caused by the diamagnetic parallel heat flow in 
GS2. The heat flow is similar in size to the heat flow for the Cyclone Base Case in figure 2, but 
the poloidal variation is not considered for simplicity. We consider the heat flow with (a) no 
radial gradient, (b) negative gradient and (c) positive gradient. 



The input of /p)’ 9 ™ 

GS2 results 


q\\/vuPi 

(dq\\/dr){a/v ti Pi) 

rUnt Vtj 
Qi Ro 

a (Ikn± m.) 

G V Qi Ro) 

(a) Constant q« 

0.1 

0 

-0.005 

0.016 

(b) Positive (dq^/dr) 

0 

0.5 

0.057 

0.018 

(c) Negative (dq\\/dr) 

0 

-0.5 

-0.067 

0.018 


the momentum flux due to the neoclassical potential (f)\° and the distribution function 
j-hg^other are non-negligible. The different results in Table 1 and |4] can be used to explain 
why the total momentum flux is reversed as the collisionality increases, as found in 
[22]. For the Cyclone Base Case, as the collisionality increases (z/* = 0.043 —* 4.3), 
the momentum flux changes its sign from negative to positive ((II mt/Qi){vti/Ro) — 
—0.022 —> 0.049). The sign reversal is due to several combined effects. 

One is the increase of the positive momentum flux due to the increase of the 
diamagnetic particle flow. Because strong collisions reduce the temperature driven 
diamagnetic flow, whose direction is negative in this case, the total size of the 
diamagnetic flow and its radial shear increases {{yi^^Rq/vu) = 0.087 —> 0.180 and 
dfop'dRo/vt^/dr = —0.152 —» —0.279). Compare the red graphs in figure 5-(a) and 
figure 2-(a). For this reason, the positive momentum flux + IIv^ increases, 

as shown in Table [5]- (a) and (b) compared with Table 2-(a) and (b). The second 
effect is the reduced negative momentum flux due to the reduced poloidal dependence 
effect ((Hf nt /Qi)(vti/Ro) = —0.012 in Table 2-(c) —> 0.001 in Table |s]-(c)). Compare 
the blue graphs in figure 5-(b) and figure 3-(b). This reduction is also due to the 
reduced temperature driven diamagnetic flow. The third effect is the reduced negative 
momentum flux due to the heat flow Qi){v t i/R q) = —0.075 —* —0.045), which 

is due to the reduced size of the heat flow. Compare the green graphs in figure 5- 
(a) an figure 2-(a). There is also a contribution due to (ftp 1 and f^ other f or moderate 

collisionality, giving the positive momentum flux + B°*t er in Table @ 

4-2.1. Neoclassical potential. To find the contribution of the neoclassical potential, we 
calculate the time averaged momentum transport in GS2 simulations by eliminating 
the neoclassical distribution function f b9 ' d = 0 but keeping the potential from NEO 
(f >\ 9 . Table [i] shows the momentum transport in terms of various collisionalities. The 
size of the neoclassical potential varies significantly with collisionality |36i 137] . As the 
collisionality increases, the size of the neoclassical potential increases and gives a larger 
momentum transport. The momentum transport due to the neoclassical potential 
in Table [6] does not increase monotonically with increasing collisionality because the 
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(a) 



(b) 



9 [Radian] 


Figure 5. The normalized diamagnetic particle flow in terms of poloidal angle calculated by 
NEO for the Cyclone Base Case with p* = 0.01 and moderate collisionality = 4.3. (a) 
The parallel particle flow (red) V\\ td /vu and the parallel heat flow q\\/vuPi (green), (b) the 
parallel particle flow (red) V\\ td = (f^, P + ^ <p,t)I/B , the pressure-driven flow (green) f 
and the temperature-driven flow (blue) Q, V) tI/B = ( K/n)B . Here, Q v , p Ro/vu = 0.167 and 
(. K/n)B 0 /v ti ~ 0.013. 


Table 4. The radial flux of the toroidal angular momentum Hint in the non-rotating state 
(fi v = £ltp,d. + &lp,e = 0), divided by the ion heat flux Qi given by GS2 for the Cyclone Base 
Case with p* = 0.01, d(R 0 /L n )/dr = 0, d(Ro/Lx)/dr = 0 and moderate collisionality = 4.3 
for (a) full neoclassical function and potential, (b) only the particle flow piece f \ 9 ?, (c) only 
the parallel heat flow piece f \ 9 ^, (d) only the neoclassical potential and (e) only the remaining 
piece fbgMhef jj erej a ^rkoi j s the standard deviation of the normalized momentum flux 
from the time averaged value due to the turbulent fluctuations. 




Size 

GS2 results 




II int Vti 

n ( Hint VtA 




Qi Ro 

a { Qi Ro) 

(a) 

IT int 

rbg,V\\ pbg,q\\ rbg, other 

J 1,2 ' Jl,i Jl,i 

0.049 

0.035 

(b) 

n v " 

V\\, d in figure 5-(a) 

0.058 

0.027 

(c) 

n 911 

LL mt 

q\\ in figure 5-(a) 

-0.046 

0.033 

(d) 

n^ 9 

LL mt 

e\/{{<t> b i 9 ) 2 U/Ti^ 0.009 

-0.011 

0.036 

(e) 

yrother 
11 int 

other than Vj| ^ and q\\ 

0.034 

0.028 
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Table 5. The normalized intrinsic momentum flux caused by the diamagnetic particle flow 
in GS2 for the Cyclone Base Case with p* = 0.01, d(Ro/L n )/dr = 0, d(Ro/LT)/dr = 0 
and moderate collisionality v* = 4.3. Here, IlJ/ t has three contributions from (a) the different 
momentum pinches , (b) the different momentum diffusivities and (c) the different 

dependence on the poloidal angle Hf nt . 


Characteristics of non-Maxwellian f^ d GS2 results 



^ip,dRo 

RqCL 

9) 

n int vu 

a (cu s tv li \ 

G V Qi R o ) 


Vti 

dr vti 

Qi Ro 

(a) ntf 

0.180 

0.0 

0.0 

-0.024 

0.022 

(b) 

0.0 

-0.279 

0.0 

0.066 

0.030 

(c) nL 

0.0 

0.0 

Fiq. 5-(b) blue graph 

0.001 

0.029 


d> ba 

Table 6. The radial turbulent flux n^ t in GS2 due to the neoclassical potential in the cyclone 
case for various collisionalities v*. The NEO code is used to find the normalized neoclassical 
potentials eyj{((/i 9 ) 2 )ip/Ti at r/a = 0.54. For the results in this table, we assume d(Ro/L n )/dr = 

0 and d(Ro/Lr)/dr = 0. The radial gradient of the potential ( ae/Ti)(d\J ((bi S ) 2 )p/<9r) is 
obtained from NEO results at two different radius at r/a = 0.53 and r/a = 0.55 using finite 
difference. 

NEO results GS2 results 

^ ey/i^h/Ti {ae/md^t^YU/dr) 


(a) 0.043 

0.0005 

0.0001 

0.002 

0.019 

(b) 0.12 

0.0012 

0.0000 

0.003 

0.020 

(c) 0.43 

0.0027 

-0.0006 

-0.011 

0.026 

(d) 1.2 

0.0048 

-0.0008 

-0.015 

0.033 

(e) 4.3 

0.0092 

0.0070 

-0.011 

0.036 


n, 


'.t Vti 


Qi Ro 


neoclassical potential at different radial and poloidal position has a different dependence 
on the collisionality. 

The effects of neoclassical potential (j)\ 9 (ij ;, 6) on the momentum transport are 
determined by two terms in the higher order gyrokinetic equation G3- One term 
gives the E x B drift of the turbulent fluctuations, —(c/B)\7(j) J / 1 x b • V// 6 . The other 
term corresponds to the parallel neoclassical electric field changing the energy of the 
particle, 


Zi,e 


bg 


-U||b ■ 

rrii 


dft 


tb 


Z,e 


Dyb • V 0 


W dfl 


tb 


(48) 


d£ rn t 09 d£ 

The energy change in the neoclassical potential due to the magnetic drift can be 
neglected because it is smaller than other terms 
vm • Vc/q 9 <C f|jb ■ V </\ g ). 


(i.e. |v M • < |v M ■ V0o 9 | and 


4-2.2. Other neoclassical contributions. To calculate the neoclassical contributions that 
are not the neoclassical particle and heat flow and the potential, we use GS2 simulations 
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Table 7. The radial turbulent flux n°*(f er in GS2 due to f b ^ other f or various collisionalities 


in the same cyclone case as Table [6] 





GS2 results 

R* 


n i 7 i Er 

vu a ( n in?'’ r Vti \ 


Qi 

Ro ° \ Qi Ro) 

(a) 

0.043 

0.001 

0.016 

(b) 

0.12 

-0.004 

0.018 

(c) 

0.43 

0.005 

0.021 

(d) 

1.2 

-0.006 

0.028 

(e) 

4.3 

0.034 

0.028 


= 0 

and / 

bq,other 1 • l 

yi , which IS 


function from NEO once the particle flow and the parallel heat flow pieces are subtracted, 


M,V\\ 


M,q\\ 


fi 9 i ~ fu V " ~ f'T" ■ Here, / 1 &9 / V|1 and /j^’ 911 are reconstructed from equations (|25|) and 
(33) using V\\ yd and gy of from NEO. 


As the size of f h ^ other increases with collisionality, the contribution II"^ er generally 
increases with collisionality. Table 7 shows the momentum transport in terms of various 
collisionalities. The contribution W^ er is approximately zero for low to moderate 
collisionalities. At the highest collisionality considered, II °fy er takes on an appreciable 
positive value that impacts the overall momentum flux. This contribution is related to 
n£ t because f^ other contains the poloidallly-dependent density and temperature that 

drive (f)\ 9 . As shown in section 
collisionality. 


4.2.1 


and Table 6, the contribution IIA increases with 


ti 3 


5. Discussion 

In this paper, the turbulent momentum transport in the absence of total toroidal flow 
is investigated. We have studied contributions to the momentum transport due to 
neoclassical effects. The momentum transport has a piece due to the difference in 
the momentum pinches IIand a piece due to different diffusivities Ilf^ between 
the diamagnetic particle flow and E x B flow. Even for diamagnetic particle flow 
that does not have flux surface averaged toroidal angular momentum, the different 
poloidal dependence between the temperature driven diamagnetic flow and the pressure 
driven diamagnetic flow can result in the momentum transport Tl 0 int . Also, there 
are contributions due to the neoclassical parallel heat flow II^ t and the neoclassical 

J>9 

potential II. Those contributions are evaluated for an example case using higher 
order gyrokinetics corrected to retain the neoclassical flow effects. We have shown that 
the relative importance of all these contributions changes with collisionality. 

The size and the sign of each contribution evaluated in section [4] for the Cyclone 
Base Case using a set of plasma parameters may not be generalizable to the whole 
parameter phase space. For the numerical results in this paper, d(R 0 /L n )/dr = 0 and 
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d{R 0 /L T )/dr = 0 are chosen, however changes in these second derivatives can modify the 
contributions to the momentum transport significantly due to the changes in dVna/dr 
and dq\\/dr. Nevertheless, clasifying the components of momentum transport by their 
driving mechanisms as is done in this paper is useful to analyze and interpret momentum 
transport with other plasma parameters. 
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